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Abstract 

In two previous papers two evolution equations for the vortex line 
density L, proposed by Vinen, were generalized to rotating superfluid 
turbulence and compared with each other. Here, the already general- 
ized alternative Vinen equation is extended to the case in which coun- 
terflow and rotation are not collinear. Then, the obtained equation 
is considered from the viewpoint of non-equilibrium thermodynamics. 
According with this formalism, the compatibility between this evo- 
lution equation for L and that one for the velocity of the superfluid 
component is studied. The compatibility condition requires the pres- 
ence of a new term dependent on the anisotropy of the tangle, which 
indicates how the friction force depends on the rotation rate. 

1 Introduction 

Quantum turbulence is described as a chaotic motion of quantized vortices 
in a disordered tangle, which in some occasions shows diffusive or undulatory 
behavior pQ-[7]. The measurements of vortex lines are described in terms of 
a macroscopic average of the vortex line length per unit volume L (briefly 
called vortex line density and which has dimensions length^ 2 ). 

It is well-known in literature that when one apply a rotation to a sample 
filled of helium II, a vortex array aligned to the rotation axis is formed. In 
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the steady state L tends to 20,/ k, fi = |0| being the modulus of the angular 
velocity and k = h jm the quantum of vorticity (m the mass of the 4 He atom 
and h Planck's constant, k ~ 9.97 10~ cm 2 /s). 

In the experiments on thermal counterflow, instead, an almost isotropic 
tangle of quantized vortices is created, with L 1 / 2 = jV/k, 7 being a di- 
mensionless coefficient which depends on temperature and V the modulus 
of the counterflow velocity defined below. These experiments are performed 
by means of channel filled of helium II and heated at one end. Using the 
two-fluid model, since only the normal fluid component carries entropy and 
heat flow, its motion will be from the heated end to the opposite end. At the 
same time, to conserve the mass the superfluid component must counterflow 
towards the heater. Thus, a relative counterflow between the normal fluid 
and superfluid components is established, characterized by V = V n — V s 
(V n and V s being the velocities of the normal and superfluid components 
averaged in a small volume A). When counterflow velocity is higher than a 
critical value then quantized vortices will be created. 

Neglecting the influence of the walls, an evolution equation for L under 
constant values of the counterflow velocity V was formulated by Vinen [8] 

^ = aVL^ 2 -f3KL\ (1) 

a and f3 being dimensionless parameters. The steady state solutions of 
equation ([1]), apart from L = 0, is L = (aV/ fin) 2 . 

There exists another version of (PQ) , the so-called alternative Vinen equa- 
tion [9]-[lQ] 

^=A^L-^L 2 , (2) 
at k 

which is also admissible on dimensional grounds. The steady state solution 
is L = A±V 2 1 fin 2 , in agreement with the experimental results in completely 
developed turbulent regime. 

The natural question which arises is what happens when the combined 
situation of rotation and counterflow exists. Supported also by experimental 
results [HI [12], this has addressed many recent studies [13]-[14] in this 
direction, some of which focusing the attention to extend Vinen's ideas to a 
wider range of situations [15|— 118|. 

One of the open problems on this arguments regards the behaviour of 
the turbulence when the rotation axis and the counterflow direction are 
neither parallel nor orthogonal, two distinguished cases already taken into 
account experimentally \U\ I12j. The parallel case, considered by Swanson 
et al. in was already investigated in Ref. [15| - - where the authors 
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extended Vinen equation (pQ) to this more complicated situation — and in 
Ref. [18], where instead the alternative Vinen equation ([2]) was extended. A 
comparison of the theoretical studies and the experimental results was only 
admissible for the experiments performed by Swanson et al. A motivation 
for the analysis of the present paper could be the behavior of a sphere filled 
of superfluid, which is submitted both to a rotation around its own axis (such 
as neutron star) and to a radial heat flux. Thus, the interaction between 
rotation and heat flow forming an arbitrary mutual angle arises in a natural 
way in these systems (see figures [2]) . 

Because of the lack of experiments on rotating counterflow situation 
with V and f2 forming a generic angle, some theoretical studies were made 
in Refs. [161 Ej extending the modified Vinen equation, proposed in [15], to 
this case. There, a thermodynamic analysis to determine possible coupling 
terms between the evolution equations of L and V was performed. 

In Ref. |18| , the extension of Vinen equation ([2]) to the case of V parallel 
to ft was considered in order to explore whether these further arguments 
are enough to establish which of both equations, (P) or ([2]), is more suitable 
to describe actual experimental results. Of course, to have a satisfactory 
version of the macroscopic alternative Vinen equation (the same arguments 
are also valid for the extension of Vinen equation), the respective coefficients 
of all terms will have to be microscopically calculated and found to coin- 
cide with macroscopic observations. These perspectives are still far from 
the present abilities because of the enormous difficulties arising to model a 
rotating tangle. Thus, a combined effort in macroscopic and microscopic 
perspectives seems a reasonable and promising way to proceed. 

The main aims of this paper is to extend the modified alternative Vi- 
nen equation, already proposed in [18J, to include situations where the angle 
between V and fi is not assigned (Section[2j). Furthermore, the compati- 
bility between the generalization of the alternative Vinen equation and the 
evolution equation for the velocity V s is studied from the non-equilibrium 
thermodynamics viewpoint (Section[3]). The new equations (equations ([6]) 
and (|23p ) model a class of physical phenomena where very few experiments 
are available. However, the model appears to be confirmed by its ability to 
describe related physical phenomena as shown in the following sections. It 
is to note that the results of the present paper need the support of the ex- 
periments to be confirmed and improved. In particular, the model could be 
applied to the neutron stars which are rotating compact self-gravitating ob- 
jects whose bulk consists of a superfluid neutron liquid at the huge density, 
with a very small concentration of protons and electrons. 
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2 Contribution of rotation to the evolution equa- 
tion for the vortex line density: general case 



In Ref. [II], Swanson et al. and, in Ref. [12], Yarmchuck and Glaberson ap- 
plied a heat flux to a rotating sample showing a complex interaction between 
processes of formation and destruction of vortices. Particularly, in the exper- 
iment performed by Swanson et al. (fi and V collinear), they observed that 
the effects of V and are not additive: in fact, for sufficiently high values 
of fl, when the influence of the walls can be neglected, the total vortex line 
density is lower than Lr + Ljj, Lr and Ljj being the values of L in steady 
rotation and in steady counterflow superfluid turbulence, respectively, given 
by 

L = L R = — , £^ = 7—, (3) 

K K 

and the deviation increases with V and f2. Therefore, the rotation facilitates 
the vortex formation, in the absence or for small counterflow velocities, 
but it hinders their lengthening for high values of V and f2. Furthermore, 
they observed that there are two critical counterflow-rotation velocities V c \ 
and V C 2, proportional to il 1 / 2 (V c \ = C\\/Tl, V C 2 = C^Vfi, with C\ = 
0.053cm sec" 1/2 , C 2 = 0.118cm sec~ 1/2 ). For V < V c2 the length L per 
unit volume of the vortex lines is independent on V and agrees with the 
expression (JSH) , but with a slightly different proportionality constant in the 
range V& < V < V C 2- Finally, for V > V C 2, L increases and becomes 
proportional to V 2 at high values of V. 

The experiments performed, instead, by Yarmchuck and Glaberson in 
[12j refer to an array of vortices, caused by a rotation Q, which suffer the 
influence of an orthogonal counterflow. They measure the gradients of tem- 
peratures and chemical potential which confirms the presence of two critical 
velocities for the counterflow velocity, the first corresponding to a depinning 
of the vortices from the walls and the second one corresponding to the tran- 
sition to the turbulent state. The critical velocities are proportional to vfi 
as pointed out by Swanson et al. in their experiments. 

As mentioned in the Introduction, in the regime of high rotation (0.2 
Hz < VL/2-k < 1.0 Hz and < V 2 < 0.2 cm 2 /s 2 ) and fi parallel to V, 
equation ([T]) has been generalized into [15] 

dL n t2 

— = —pnL + a.\ 



rl/2 

L 1 — mi — ;= 



VL + fo 



rl/2 

L 1 — ni2 — — 



'kUL, (4) 
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where m\ = fi^/oti and m<i = Pi/fa, and j3, a\, f3\, fa and /J4 are coefficients 
depending on the polarization of the tangle, which was supposed to be a 
function of and V. As observed in [15], equation Q describes some of the 
most relevant effects observed in the experiments of Ref. [11] . 

Since experimental results in pure counterflow are also described by the 
alternative Vinen equation, in [18] a generalization of it to rotating counter- 
flow was proposed following the main lines of Ref. |15j . This generalization 
was carried out assuming that the counterflow velocity and the rotating axis 
were parallel one to each other in order to compare theoretical results with 
experiments by Swanson et al. This equation assumed the following form 



— = -$kL 2 + A x 
at 



r jr 


V 2 




L — V\ — 


— + B, 


L - v 2 — 


K 


K 


K 



n. 



(5) 



and it is valid only for V collinear to fl. In Eq. ([5]) coefficients j3, A±, u\, B\ 
and U2 were supposed to depend on the polarization of the tangle, i.e. on 
fl and V. 

Now, following the general line of the paper [17] equation ([5|) is written 
to include the general situation in which vectors ft and V are not necessarily 
collinear. Looking at the terms of equation ([5]), the only one which needs 
experimental examination is that proportional to £IV 2 , because in the gen- 
eral case it could depend not only on the absolute values of the two vectors 
ft and V but also on the angle between them. A mathematically consistent 
version of ([5]) including both alternatives is 



dL 
~dt 



(3kL 2 + L 
AiVvx 



A, 



V • U • V + Bift ■ U • ft 



ft 



ai Vfi + a 2 fiv) ■ V + ^^-ft U ft 



(6) 



where U is the second order unit tensor, at and a% are two coefficients 
such that a\ + ai = 1, Vfi and f2V are diadic products between V and ft 
such that the first depends only on the absolute values of two vectors V, fl, 
whereas the second one depends also on their angle, cos 2 (Vfi). 

As pointed out above, equation ([6]) includes different situations and it 
becomes equation ([5]) when one considers a rotating cylinder heated from 
the lower basis or the higher basis. A physically interesting situation arises 
when two vectors are not parallel, for instance when one consider a rotating 
cylinder heated from the wall, which means that V and ft are perpendicular 
as in experiment performed in Ref. [T2j (these experiments will be discussed 
in the Conclusions). In the general case, it is important to establish which 
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term, a± and has to be included into equation ([6]). Since up to now we 
are not able to answer to this question, both coefficients a\ and a 2 have to 
be considered. 

The stationary solutions of equation ([6]) are obtained requiring dL/dt = 
0, that is 



-f3nL 2 +L ( — V 2 + B t Q 



2 ( a i + °2 cos (" v ) J H — " 

Kj \ /Hi 

The case fi parallel to V, that is the stationary solutions of equation jZJ) 
with cos 2 (f2V) = 1 and a\ + 02 = 1 was investigated in [18]. The main 
results of the paper [18] are reported in order to extend them to the general 
case dealt with in this paper. 

In the steady state (L and V constant), under the constraint 

v 2 = v x -v\^, (8) 

the stationary solutions of equation © can be written 

L = l\=u x -, (9) 

II Ai V 2 (Bi \ 

Solutions ([9]) and (|10p are two families of straight lines in the plane (V 2 , L), 
as plotted in Fig. [TJ the first of them (equation (J9|)) parallel to the V 2 axis 
and the second one (equation (flOl) ) with a slope independent of fi. Applying 
the stability analysis to these solutions, solution ([9]) is stable if V is lower 
than 

K2 = 1 fire, 11 



0. 



while, for V higher than V c i the solution (|10p is stable. The obtained 
value of the critical velocity (|11|) corresponds to the interception of the two 
straight lines (|9|) and ()10p in Fig{TJ It also corresponds to the second critical 
counterflow-rotation velocity observed in the experiments of Ref. |llj . As 
we see, this critical velocity scales as \/0, in agreement with experimental 
observations. 

To describe the existence of the first critical velocity V c \ in the experi- 
ments, and also the small step in L, in Ref. [18] v\ was assumed dependent 
on $7 and V as 



Vl = A \ 1 - B tanh 



(12) 
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There, A = 2.018, B = 0.0089 and C = 0.355 were found in such a way 
that for V 2 <C V 2 X = lyktt, it results v\ ~ A - BA = 2 and for V 2 > V^, 
v\ = A + BA = i/ max = 2.036, while the constant C is related to the critical 
value Vci of the counterflow velocity by V c 2 = ^kQ, whose experimental 
value is V c i = 0.053\/^ cm sec -1 / 2 . In (|12p . N' is a phenomenological 
coefficient characterizing the rate of growth of L near V& and for it the 
value N' = 22 proposed by Tsubota et al. in Ref. [2] was assumed. 

Now, let us investigate on the stationary solutions of equation ([6]). Thus, 
assuming that relation (jSJ) is still valid when one substitutes v\ with v\ = 
v\(ax + a 2 cos 2 (f2V)), that is 



v 2 = v\ - (13) 
B\ 



then (UJ) has the following stationary solutions 

L = L\ = D\— = ia\ + a 2 cos (12V))— , (14) 

fC K 

A 1 V 2 : (B 1 _\n A 1 V 2 ,(B 1 , , 2/ft V ))^ n 

(15) 

The stability analysis furnishes the same conclusions of the parallel case, 
namely L 2 is stable for V > V c2 , V c2 being the second critical velocity defined 
by 

V r n = — ilK = ilK (16) 

A 1 vx-v2 A x 

and depending on the angle between Cl and V through v\. Looking at the 
first solution L\ one notes that 

ai < (ai + 0,2 cos 2 (fiV)) < (oi + a 2 ) = 1 

which means 

v\ being the coefficient of the stationary solution L% in (j^J) and plotted in 
Fig.[TJ An equivalent relation could also be written for L 2 , but, differently 
to L\, coefficients A\ and B\ could depend on the anisotropy of the tangle, 
namely on the angle between V and ft. 

In Ref. [15] . the dimensionless quantities appearing in the evolution equa- 
tion ([5]) are determined by comparison with the experimental data in the 
steady state, obtaining 

A B 

-j = 0.0125, -j = 3.90, z/i = 2.036, u 2 = 0.97. (17) 
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Instead, in the case of pure counterflow A\j '/3 = 0.156, which confirms the 
dependence of the coefficients on the anisotropy of the tangle [18] . By (fT7|) . 
solutions (|9|) and (fTUj) become 

L\ = 2.036-, and 4 = 0.0125^- + 1.86-. (18) 

K K K 

These solutions are plotted in Fig.[TJ where the agrement with the experi- 
mental data is shown. 




0.05 0.1 0.15 0.2 



Figure 1: Comparison of the solutions (118p (continuous line) with the ex- 
perimental data by Swanson et al. [UJ. Figure from Ref. [18] . 

An interesting problem is to establish which equation, either the one 
based on the usual Vinen equation or the other one based on the alternative 
Vinen equation, fits better the experimental data obtained by Swanson, 
Barenghi and Donnelly [TT]. This comparison was made in Ref. [18] : the two 
stationary solutions in the range V < V C 2 represent the same straight line 
in the plane (L, V 2 ), whereas in the range V 2 > V 2 2 it is not the case. From 
the direct comparison with the experimental data — by finding the errors 
between the stationary solutions of the two models and the corresponding 
experimental values — in Ref. [18] the following statement was established: 
the stationary solution of the alternative Vinen equation approaches better 
the experimental data (for V 2 > V 2 2 ) than that of the usual Vinen equation. 

In the same paper the unsteady behavior of the solutions of these equa- 
tions was also explored. Both equations exhibit remarkable differences, such 
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as the solutions of the alternative Vinen tend much faster to their steady- 
state values, and this difference depends on the value of the counternow 
velocity. Even though detailed experimental data on this unsteady behavior 
lacks, according to Swanson et at the time required to reach the steady state 
is less than 10 minutes — when the counternow velocity V is slightly above 
the critical velocity V C 2 and it increases between two consecutive experimen- 
tal values (see pag. 191, Ref. [II])- According to the results of Ref. [IS] , 
the temporal scale of the solution of the usual Vinen equation was closer to 
the observations than the temporal scale corresponding to the alternative 
equation, which tends too fast to the final result. Thus, it seems that the 
usual equation is preferable on these grounds. 



3 Thermodynamic analysis 

In this Section the general lines of Refs. |10 |, I16 |, ITT] are followed to determine 
an evolution equation for V s which is consistent with equation ([6]). This 
situation is of special interest when V s is not kept constant, as it is usually 
assumed in ([I]) or ([2]) — or in ([4]) and ([5]). This more general model may 
be useful to study the dynamic behavior of the vortex tangle, when the 
counternow velocity V varies in a greater range than that considered in 
Ref. [18]. According to the formalism of nonequilibrium thermodynamics, 
evolution equations for V s and L can be obtained by writing <JV s /dt and 
dL/dt in terms of their conjugate forces — /? S V and ey. This is motivated also 
by the paper of Nemirowskii et al. |10| . where for the entropy density s of 
the superfluid in the presence of vortex lines they considered the differential 
form 

ds „ T gTV, dL 



where 



p s V 



dt 



du 



9V S 



-PsV 



dt 

du 
dL 



ev Tt' 



Ps^ 

4tt 



In 



(19) 



(20) 



being u the internal energy density and ey the contribution to the internal 
energy per unit length of the vortex line (ao is the dimension of the vortex 
core, which is very small, of the order of one A) [TJ. 

First, recalling that in the presence of pure rotation (which produces 
an ordered array of vortex lines parallel to the rotation axis) the evolution 
equation for V s is [7J ETJ CE] 



Pn B 



+ 2Q x V s + i = —L—— !1xflxV + — ft x V 

dt pax B 



B' 



(21) 
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where a = 2/k, iq is the inertial force, B and B' are the Hall-Vinen dimen- 
sionless coefficients describing the interaction between the normal fluid and 
the vortex lines |19] . 

In the presence of pure counterflow, the evolution equation for V s can 
be written [TU] 

^JL = El^ALV, (22) 
dt p ' v ' 

where A = 2B/?>a [7J. Since in the steady state L is proportional to V 2 , 
the friction force on V s (which is opposite to V), will be proportional to 
V s , a result known as Gorter-Mellink law [7j [2j [10]. Note, however, that 
this V 3 dependence is valid only at the steady state, whereas ([22]) may be 
also applied in unsteady situations, in absence of rotation. Here we will be 
interested in the form of the friction force in unsteady situations including 
rotation and, in more general terms, in the general form of the evolution 
equation for V s . 

In the simultaneous presence of counterflow and rotation, terms as those 
in the right-hand side of (I21D and (j22[) must be included in the evolution 
equation for V s , as well as an additional contribution F coup i = F C (V, ft) 
due to couplings between counterflow, rotation and superfluid velocity. The 
coupling term is found similarly to the Onsager's formalism of non-equilibrium 
thermodynamics p~Ql [EU [17] 

^± + 2Qx V s + i = ^L?n.V + F c (V,ft), (23) 
dt pa 

where tensor 

2 / \ B' 

n = (i-6)-u + 6(u-nn) +b'— w n (24) 

3 V / B 

was introduced and studied in Ref. [20]. In Refs. [16] 117] . for the sake of 
simplicity, b = b' was chosen. Here, b and b' are numerical parameters related 
to the anisotropy and polarization of vortex tangle, describing the relative 
weight of the array of vortex lines parallel to CI and the isotropic tangle: 
b = b' = means an isotropic tangle corresponding to pure counterflow 
and b = b' = 1 means an ordered array corresponding to pure rotation. 
The microscopic physical interpretation of coefficients b and b' was made in 
Ref. |20| . In particular b is linked to the anisotropy of the tangle and b' to 
the polarity. 

Our aim here is to include equations © and (j23|) into a common ther- 
modynamic framework, and study possible couplings between them. But, 
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because of the origin of the term QV 2 in equation ([5]) is unknown, as pointed 
out in Section^ here special emphasis is given to the two opposite cases: 
a% = 1 and a-i = 0, and a\ = and a>2 = 1. Of course, an intermediate 
situation where both coefficients are non zero is not excluded and could be 
an interesting physical situation which will not be considered for the sake 
of simplicity. Future experiments could urge us to spend further efforts on 
these intermediate cases. 



3.1 QV 2 interpreted as Vtl VCl V (ai = 1 and a 2 = 0) 

In this subsection, term QV 2 is assumed to come from Vfl • Vfi • V which 
depends only on the absolute values of the two fields. In this spirit, equa- 
tion ([6]) for L with a\ = 1 and 02 = is written in the second line of 
the following system (I25D . The second part of the evolution equation for 
V s is built up by means of the Onsager-Casimir reciprocity relation, as in 
Refs. [MCE]. The result is 



[DV S - 




dt 




dL 




- dt - 





1_A 

Ps 



PPs a 



±Al\ L 

Ps K L 



M^-*g]V -±[P*I?-Bi{L-vfi)ti\ 



where 



DV, 



+ 2fl x V s + i . 



ev ' 
(25) 

(26) 



dt dt 

Therefore, in the presence of counterflow and rotation the equation for 
dV s /dt is 



dV 



, +2fixV s + io = L^^n.V+^^ 

dt pa p s k 



ev A! 



L — V\ — 
K 



V. 



(27) 



From a comparison with (|23D the following expression for the coupling fric- 
tion force is obtained 



F — 

r coupl 



ey_ Ai 

Ps K 



L 



n 



V. 



(28) 



This term is proportional to V, and is therefore a friction term (recall that 
V is opposite to V s ). Note the presence of a negative term proportional 
to and independent on L, which means that the rotation of the vortex 
tangle reduces the friction. This may be interpreted as a consequence of 
the alignment of the vortices to the rotation axis due to ft, in such a way 
that they have less resistance to the flow (recall that the friction is due to 
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vortices orthogonal to the counterflow velocity, but not to those parallel to 
it). 

Substituting (|14p and (|15p in the off-diagonal term in the matrix in (|25p , 
one obtains the following expression for the coupling force 



tt — 

r coupl — 



eyAi 

Ps K 



coupl 



Ps k 



Li 



n 



V 



v = o, 

Ps ^ 



for V < V, 



c2, 



(29) 



K 2 2 )V, 



for V > V c2 . (30) 



As a consequence, in a steady situation (and also in almost steady state) for 
V < V C 2 the coupling force is absent (as in pure rotation) while, for V > V c2 
— when the array of rectilinear vortex lines becomes a disordered tangle — 
the additional term (|30p proportional to V appears. Thus, V c2 indicates the 
threshold not only of the vortex line dynamics but also of the friction acting 
on the velocity V s itself. 



3.2 QV 2 interpreted asVfi-fiV-V {a x = and a 2 = 1) 

Here the same lines of the previous subsection are applied to the opposite 
case a± = and a 2 = 1, which means that term QV 2 comes from Vfi-fiV-V, 
depending also on the angle between the two vectors. So, writing again 
equation ([6]) with a± = and a 2 = 1 in the second line of the following 
system (i31~|) , the second part of the evolution equation for V s is built up by 
means of the Onsager-Casimir reciprocity relation 



DVs 
dt 
dL 
dt 



Pa K 



pps a 

L-L 



Ps K 



L-L-i 



V 



V -A-[p K L 2 - Bl {L-u 2 %)n] 



-PsV 



(31) 

where L\ = v\— cos 2 (f2V). In this case, instead, the equation for dV s /dt in 
the presence of counterflow and rotation is 



—— + 2fi x V s + 1 = L ^11 ■ V H 

at pa p s k 



L-z^i-cos 2 (fiV) 



V. (32) 



From a comparison with (|23p . the following expression for the coupling fric- 
tion force is obtained 



coupl 



eyAi 

Ps K 



L - ^i-cos 2 (fiV) 

K 



V. 



(33) 
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Note that it is equal to (j28|) unless cos 2 (fiV). So, all the observations below 
relation (|28p are still valid. Note also that the negative term, proportional to 
Cl and independent on L, depends explicitly on the angle between V and Cl 
and it is maximum for V parallel to Cl whereas it is zero for V orthogonal to 
Cl. The friction term has always the same direction of counterflow velocity 
V, and in the situation V_Lfi depends only on V and L but not on the 
angular velocity Cl. 

By using expressions (fTl")) and (fT5j) , the coupling term ([33]) in the steady 
state becomes 



coupl 



coupl 



_ evAi 

Ps « 


Li 


ci 

- v x — 

K 


_ evAi 

Ps K 


L 2 


Cl 

- v\— 

K 



v = o, 



V 



for V<V c2 , 
(34) 

3 (V 2 -V£)V, iovV>V c2 . 

(35) 



Here, the conclusions are the same pointed out below relations (|29j) and (|30j) . 
and, therefore, also in this case V^2 indicates the threshold of the friction 
force acting on the velocity V s itself. 



4 Conclusions 

The possibility of at least two reasonable evolution equations for the vortex 
line density L in superfluid turbulence was known since the early days in 
which Vinen proposed them. This equation was generalized to the combined 
situation of counterflow and rotation, with V parallel to CI, into (jl]) and 
([5]), which was compared with the experimental data given by Swanson et 
al. . In [T7] Vinen equation (jlj) was extended to include the case V and CI 
with different directions. Here, following the same ideas of that paper, an 
extension of the modified alternative Vinen equation (J5]) has been proposed 
in©. So, this equation is valid when V and Cl have indefinite directions, 
becoming the evolution equation (JEJ when these vectors are parallel. The 
stationary solutions of this extended equation differ from that of parallel case 
for the coefficient v instead of v. So, from the theoretical point of view one 
can establish that an higher angle between Cl and V implies a lower value 
for the steady state L%, unless future experiments will prove that v = v. 

The case V orthogonal to Cl was experimentally carried out by Yarm- 
chuck and Glaberson in Ref.[12j. They arranged a pair of horizontal parallel 
glass plates forming a closed channel of rectangular cross-section closed at 
one end with a heater nearby, and open at the other end to the liquid helium 
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bath. The channel was of large aspect ratio: the height being 0.5 mm, the 
width 1.4 cm, and the length 5.5 cm. The channel is rotated around its 
vertical axis, orthogonal to the direction of the heat flux, in such a way that 
the counterflow velocity V is orthogonal to angular velocity ft. The authors 
check temperature and chemical gradients as functions of heater power and 
rotation speed. The results were linear regimes, in which the temperature 
and chemical gradients increase when rotating velocity is increased, and a 
critical heater power q C 2 (corresponding a critical counterflow velocity V^), 
which increase as the rotation speed increases, becoming proportional to 

Vn. 

These results points out that coefficient a\ in equation ([6|) is not zero. 
For, in the situation of the above experiment V is orthogonal to ft, namely 
the term proportional to 02 is zero because cos 2 (Vfi) = 0. On the other 
side, experimental results confirm the presence of the term proportional to 
V 2 Q. But, to establish whether the coefficient 02 is equal to zero, further 
experiments are required. 

As mentioned in the Introduction of this paper, an interesting appli- 
cation of these arguments could be a rotating spherical container (like a 
star) filled by superfluid helium which suffers the presence of radial heat 
flow (see figured]). Since irradiation comes from the center of the star, the 
heat flow will be proportional to r -2 , r being the distance from the center 
of the star. Because of the proportionality between heat flow and counter- 
flow velocity, this dependence on r requires the existence of a critical radius 
r2 (corresponding to the second critical velocity V c ?) which marks the tur- 
bulent region (r < r<z). Of course, the critical radius T2 depends on the 
angular velocity f2 from the relation (116|) . An interesting question which 
arises could be whether T2 depends on the angle 9 between V and fi, or 
similarly whether 02 is zero or not. Because of the lack of experimental 
evidence any conclusion on the presence of the coefficient a-i is not reached. 
However, the case 02 = would require a critical radius free from 6 (see 
relation (]16p ) and a spherical turbulent region of radius T2 as showed in 
figure[2k. The case (12 / 0, instead, would require a radius r2 depending 
on 6 caused by the fact that the critical counterflow velocity V C 2 increases 
when one moves from the equator to the poles of the sphere. This means 
that T2{9) has a maximum at the equator and minimum at the poles, that 
is the turbulent region is contained in an ellipsoid (see figured). Note that 
in both figures[2] the region for r > T2 is filled of vortex array which suffers 
the perturbation coming from the inner turbulent region. 

Because of the direct applications of these arguments on neutron stars, 
some other remark could be pointed out. It is to note that different forms of 
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the superfluid turbulent core have implication on the moment of inertia of 
the star which may influence the acceleration periods. These acceleration pe- 
riods are a real topic which has deserved many interest for some researchers, 
but nothing has still been proved because the interior of the neutron stars 
are not directly observable. But, some more information on neutron stars 
can be obtained because Fig. [2] can be experimentally achieved: take a ro- 
tating sphere full of superfluid, which has a small heated sphere in its centre, 
then the details of the frontier between the vortex tangle and the array of 
parallel vortices (or of a highly polarized tangle) can be measured. 

The procedure pointed out in Refs. [IT] was generalized to derive an 
evolution equation for the superfluid velocity V s , thermodynamically con- 
sistent with the proposed evolution equation for L. The form ([28]) and (f33|) 
point out — in analogy with the conclusions reached in Ref. [17J but in a 
more direct way — that a tangle of a given L and Q has less resistance to 
the flow than a tangle with the same value of L but in absence of rotation. 
This fact may be attributed to the orientational influence of the rotation on 
the vortex lines. 

In Ref. p2], the authors obtained the coupling force F conp /, between V 
and CI, in the evolution equation for superfluid velocity coupled to the gener- 
alized Vinen equation. As pointed out along this paper, they distinguished 
essentially two cases, both corresponding to ours. Regarding the case of 
subsection !3.2l — in which the investigated term of the generalized Vinen 
equation (and alternative Vinen equation) depends on the angle between V 
and fi — , in [17] the authors find F coup i ~ — L cos(Vf2)f2/\/0 which differs 
from the corresponding (|33[) obtained in this paper for the dependence on V. 
Note also that the sign of coupling force ([33]) does not depend on the angle 
between V and ft whereas that found by Jou and Mongiovf does. The cou- 
pling force (|33p is, instead, similar to the last two terms of equation (4.12) 
of paper [17] . 

As regards the case of subsection l3.ll — in which the investigated term 
does not depend on the angle between V and fi — it corresponds to case B 
of the paper [TTJ, where the authors found F coup i ~ — LVOV which differs 
from that proposed here. Again, the coupling force (j28[) is similar to the last 
two terms of equation (4.17) of paper [TTJ. 
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Figure 2: Sketch of a rotating sphere with an internal turbulent region and 
an external laminar region. The boundaries of the turbulent regions will 
be or not spherical depending on whether the critical heat flux (|16j) does 
not depend on the relative angle between heat flux and angular velocity or 
depends on it. The picture of the external region (r > r%) is only indicative, 
but a strongly polarized turbulent region could be more realistic. 
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